Frustration-induced valence bond crystal and its melting in MosSby 
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121/123 gj^ nuclear quadrupole resonance and muon spin relaxation experiments of MoaSb? revealed 
symmetry breakdown to a nonmagnetic state below the transition recently found at Ts ~ 50 K. 
The transition is characterized by a distinct lattice dynamics suggested from narrowing of nuclear 
fields. We point out that the Mo sublatice is a unique three-dimensional frustrated lattice where 
nearest-neighbor and next-nearest-neighbor antiferromagnetic interactions compete, and propose 
that tetragonal distortion to release the frustration stabilizes long-range order of spin-singlet dimers, 
i.e., valence bond crystal, which is thermally excited to the dynamic state with cubic symmetry. 

PACS numbers: 71.20.Be, 76.60.-k, 76.75.+i 
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Exotic phenomena induced by the geometric frustra- 
tion have attracted great interest, and trial to find novel 
phases in various types of frustrated lattices has been 
going on In the case of three-dimensional (3D) 

pyrochlore lattice, the symmetry breakdown to spin- 
Peierls-like states or to the valence bond crystal with 
the quenched spin degrees of freedom 0, Q is one of the 
options to release the frustration; a typical example is 
found in the case of MgTi204 On the other hand, in 
the field of correlated-elcctron superconductivity (SC), 
another phase transition prior to SC. such as the anti- 
ferromagnetic (AF) or hidden order attracts interest, 
because the particular electronic states may be a key to 
understand the unconventional SC. 

A metallic compound MosSby, with the cubic h:^Ge^ 
type structure (space group Im3m), is a superconductor 
with a critical temperature Tc ~ 2.3 K |6j]. As seen in 
the inset of Fig. [TJa), the Mo sublattice (12e) is the 3D 
network of Mo-Mo dumbbells formed by nearest neigh- 
bors (NN). Next-nearest- neighbor (NNN) bonds form an 
octahedral cage at the body-centered position. The suf- 
ficiently short NN distance (3.0 A) comparing with the 
NNN distance (4.6 A) makes us anticipate easy dimer- 
ization between NN pairs. Recently, Candolfi et al. 
regarded MoaSbr as an itinerant electron system and sug- 
gested coexistence of SC and spin fluctuations, which is a 
rare example in d electron metals. A broad maximum in 
the temperature (T) dependence of the susceptibility (x) 
with relatively large magnitude, a large electronic spe- 
cific heat coefficient, a dependence of low-T electrical 
resistivity, the Kadowaki-Wood relation between the spe- 
cific heat and the resistivity, etc. point to the presence of 
strong electron correlation. Band structure calculations 
for cubic MoaSby [1, Q, showing appreciable contribu- 
tion of Mo-4d orbitals to the density of states near the 
Fermi level, are in agreement with these observations. 
Although no phase transition other than that at Tc was 
known in fl\, Tran et al. JTd\ reported recently the pres- 
ence of another phase transition at T* ~ 50 K (in this 
article, referred to as Ts), and proposed spin gap forma- 
tion with the magnitude of 120 K associated with Mo-Mo 



dimerizaiton. Here the Mo^+ state with spin S* = 1/2 
was presumed implicitly. In fact at sufficiently high T, x 
shows a S — 1/2 Curie- Weiss behavior with a large neg- 
ative Weiss constant, suggesting strong AF interaction 
among active Mo spins. Thus the Mo electronic state is 
of interest as the system with both itinerant and localized 
characters. Furthermore MosSby is now recognized as a 
new material with another phase transition prior to SC. 
As for SC, the pairing mechanism is still controversial. 
Although most studies point to a BCS-type conventional 
SC d, M, unconventional pairing has also been 

discussed |l2l. Il3]. In it was argued that the Fermi 
surface nesting leads to SC. 

In this article, to manifest microscopic and dynamic as- 
pects of the magnetism in MoaSby, particularly to iden- 
tify the phase between Tc and Ts, we present ^^^^^^^Sh 
nuclear quadrupole resonance (NQR) and muon spin re- 
laxation (/iSR) results together with low-T x-ray diffrac- 
tion (XRD) data. We discuss possible frustration in the 
Mo sublattice by assuming AF NN and NNN interac- 
tions, and propose new exotic states realized due to the 
competition between NN and NNN interactions. 

Polycrystalline MoaSby specimen were prepared via 
reaction between liquid Sb and solid Mo. NQR spin- 
echo measurements were carried out using a conventional 
phase-coherent- type pulse spectrometer at 1.4-200 K. 
Zero-field (ZF) /j,SR measurements were made at RIKEN- 
RAL Muon Facility at the Rutherford- Appleton Labora- 
tory in the UK using a pulsed positive surface muon beam 
at 5-250 K. Low-T powder XRD was measured using a 
conventional diffractmeter (MacScience MXP) with Cu 
Ka radiation. 

Figure [TJa) shows the ^^-^/^^'^Sb NQR spectrum of 
MosSby at 200 K. Antimony has two isotopes ^^^Sb (nu- 
clear spin / = 5/2) and ^^^Sb (/ = 7/2), which yield 
two and three quadrupole transition lines, respectively, 
when the electric field gradient (EFG) is axially sym- 
metric. In cubic MoaSby, there are two nonequivalent Sb 
positions with nonaxial and axial symmetries (12c? and 
16/, respectively). We have found only five lines at 200 
K corresponding to either site. This is probably because 
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the nuclear spin relaxation rate at the unobserved site is 
enhanced due to fast fluctuations of active Mo electron 
spins. The frequency intervals between neighboring lines 
of ^^^/-'^^■^Sb are approximately equal, indicating nearly 
axial EFG at the observed site. ^^^/^^^Sb quadrupole 
frequencies are 12i/i23j,q ^ 47.7/27.2 MHz at 200 K. 

The T dependence of the ^^^Sb-2i^Q line (transitions of 
I ± 3/2) ^ I ± 5/2)) is shown in Fig. [^b). The line is 
shifted to the high frequency side and gradually broad- 
ened with approaching Tg- At ~ Ts, the line loses in- 
tensity and disappears due to highly enhanced spin-echo 
decay rate I/T2 (see the inset of Fig. [5]). Below 20 K 
the signal reappears with complicated structures (see a 
part of the data at 10 K in Fig. [Hb)). Actually ZF nu- 
clear resonance signal was continuously observed in al- 
most the whole frequency range below 110 MHz with a 
number of peaks. This indicates clearly symmetry break- 
down below Ts- The origin is either appearance of the 
internal field (due to magnetic order) , separation of crys- 
tallographic Sb sites (due to structural deformation) or 
both. Unfortunately, the spectrum at low T is too com- 
plicated and unresolved to conclude this issue because of 
line split together with broadening, as well as the pos- 
sible emergence of the signal from another Sb site. An 
early ^^^Sb Mossbauer study is not articulate either p^ . 
Hence, to determine the magnetic ground state, we em- 
ploy another microscopic probe, /iSR, which is sensitive 
to the appearance of the internal field; the result will be 
shown below. 

The nuclear spin-lattice relaxation rate 1 /Ti was mea- 
sured for both ^^^Sb- and ^^^Sb-2i/Q transitions. Nu- 
clear magnetization recovery was fitted well using the 
multi-exponential function with a sin gle parameter ex- 
pected for the case of axial EFG [T5lll6|. The ratio of 
relaxation rates for the different isotopes was estimated 
to be 123^^/121^^ ^ 3 5 at T > Tg, which is close to 
the squared gyromagnetic ratios ^'^^1^1^'^'^ In — 3-4. This 
confirms that the relaxation process is not by fluctuations 
of EFG but by those of the internal magnetic field at least 
above Tg- l/^i was measured at the peak position of the 
i2iSb-2i^Q line and plotted in the form of l/{TiT) vs. T 
in Fig. [51 The divergent behavior was observed at ~ Ts, 
although reliable values of 1/Ti could not be obtained 
near Tg due to poor signal-to-noise ratio [l7|. In fact, the 
transition at Ts is characterized by the stronger T depen- 
dence of l/Tz (the inset of Fig. H]); I/T2 ~ 10-103(l/Ti). 
Below Ts, l/(riT) decreases rapidly with decreasing T, 
being consistent with the possible existence of spin gap 
[lo| . but turns to show a TiT = constant Korringa be- 
havior below ~ 10 K, indicating the presence of residual 
density of states at the Fermi level, which is related with 
the residual x at low T. 1/Ti is enhanced just below Tc- 
Such an enhancement is known as a coherence peak and 
usually observed for the conventional BCS superconduc- 
tor with an isotropic gap [l^ . 

Now we shift to /iSR results. Typical examples of ZF- 
/iSR spectra are shown in the inset of Fig. [31 Even 
at the lowest T (5 K), no straightforward evidence of 
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FIG. 1: (color online). (a) i^Viasg^j j^qj^ spectrum of 
MosSbr at 200 K. Inset shows the Mo sublattice. (b) Tem- 
perature variation of ^^^Sb-2i/Q line. 
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FIG. 2: (color online). Temperature dependence of ^^^Sb- 
1/{T\T) for MoaSbr. Inset shows the temperature depen- 
dence of I/T2. 



magnetic order such as muon spin precession was ob- 
served. At all r, the feature of relaxation is essen- 
tially the same and shows Gaussian-like depolarization. 
The commonly used damped Kubo-Toyabe (KT) func- 
tion P^{t) = exp(-Ai)GP(A,t) with GP(A,t) = 
i -I- |(1 - A^t^) exp(-iA2t2) .^^as fit to the data, where 
A/7p is the width of the static field distribution (7^ the 
muon gyromagnetic ratio), and A is the damping rate as- 
sociated with an additional relaxation process. As seen 
below, it is likely that the KT term is dominated by the 
nuclear dipolar field and the damping part reflects the 
Mo electronic system. Figure [31 shows the T dependence 
of A. First, let us see the T-independent part in the in- 
termediate T range of ~ 70-140 K with A/7^ ~ 0.25 G. 
This tiny field is most probably due to dipolar fields com- 
ing from randomly oriented nuclear spins. Above 140 
K, A is reduced markedly. As shown in Fig. [H In A is in 
proportion to 1/T in the T range, suggesting a motional 
narrowing process caused by muon hopping among differ- 
ent sites, that is frequently observed at high T; the acti- 
vation energy is E^/kB — 910 K. With decreasing T pass- 
ing through Ts, A is enhanced strongly as seen in Fig. [31 
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FIG. 3: (color online). Temperature dependence of the ZF 
Kubo-Toyabe relaxation rate A in MosSby. Inset shows typ- 
ical examples of ZF-/iSR spectra (at 5, 55 and 150 K). Solid 
curves indicate the fit by the damped KT function. 
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FIG. 4: (color online). ZF Kubo-Toyabe relaxation rate A 
in MosSbr plotted against 1/T. Activation energies for two 
difi'erent thermal excitation processes are estimated to be 120 
and 910 K from the solid and broken lines, respectively. Inset 
shows the temperature dependence of the relaxation rate A. 



This behavior is clearly related with the transition at Tg . 
The origin of the additional internal field below Ts is the 
most essential point in this experiment. The Gaussian- 
type random distribution with small ~ 2.1 G even 
at the lowest T is hardly attributed to the electron spin 
freezing but is reasonably explained as due to the change 
in nuclear dipolar field associated with the structural de- 
formation mentioned below [l^. The T dependence of 
A at ^ Tg is again explained as a thermal excitation 
process with activation energy of iJa/^B — 120 K. This 
implies that the T variation is dominated by a dynamic 
process, giving very fruitful information on the nature of 
the transition at Tg. 

The T dependence of A is shown in the inset of Fig. ID 
The decrease above ~ 150 K would be analytical artifact 
due to the evlution of the relaxation process associated 
with the dynamic muon hopping. A has a sharp peak 
at Ts and decreases rapidly below Ts. This behavior 
is qualitatively same as that of NQR-l/Ti, and reflects 
the phase transition in the Mo electronic system. The 
quantitative difference in NQR-l/Ti and /iSR-A is most 
probably due to the difference in the coupling mecha- 
nism; the interaction between Sb nuclei and Mo electrons 
is probably dominated by the direct hyperfine coupling, 
while fields from Mo electrons to the muon site is by the 
dipolar mechanism. Different time windows of the mea- 
surements would also be appreciable. Furthermore, the 
muon stopping site may be a high symmetric position, 
say the center of the Mo cage, where the fields from Mo 
electron spins would mostly be canceled. 

The /zSR result points to the nonmagnetic ground 
state. This is in accordance with the spin-singlet dimer- 
ization suggested in TO]. To get information on the lat- 
tice symmetry below Ts, we made conventional XRD 
analyses at low T and found appreciable difference in 
the patterns below and above Ts. To see the change effi- 
ciently, we particularly show 600 diffraction peaks at 70 
and 15 K in Fig. Ella) As seen in the figure, a high-T 
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FIG. 5: (color online), (a) X-ray 600 peak of MoaSbr at 70 K 
and its split pattern at 15 K [23|. The contribution of Ka2 
radiation is subtracted by profile fitting, (b) Dispersion curves 
of the exchange integral J{q) of the Mo sublattice assuming 
negative Ji and J2. Ji = —0.05 and J2 = 0.3Ji were used 
arbitrarily. 



singlet splits into double peaks with the intensity ratio of 
'--^ 2 : 1 by crossing Ts , suggesting structural deformation 
from cubic to tetragonal (space group lA/mmm) as the 
simplest possibility; the estimated lattice parameter a is 
larger than c by '--^ 0.2% at 15 K. 

Although MosSbr is a metallic magnet, here we fo- 
cus on the local nature of the Mo electron. Taking a 
glance at the Mo sublattice, if the NN interaction Ji is 
dominant, no substantial symmetry breakdown is neces- 
sary to form diniers but isotropic shrinkage is sufficient. 
Therefore in discussing the mechanism of the transition, 
the low-T lattice symmetry is essential and NNN interac- 
tion J2 should be taken into account. Note that positive 
(ferromagnetic) J2 as suggested in [l^l would not desta- 
bilize the cubic state. Here assuming negative Ji and 
J2, we calculated eigenvalues of their Fourier transform 
J(g) for the cubic state in the scheme of the LCAO (lin- 
ear combination of atomic orbitals) approximation; see 
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Fig. IHb). Further neighbor interactions are neglected 
for simphcity. A flat dispersion of J{q) found afong the 
highest P-(D)-N branch indicates macroscopic degenera- 
cies of AF ground states 2l|; the Mo sublattice is not 
fuUy frustrated Uke pyrochfore but is apparently frus- 
trated for all AF states with q = with arbitrary ^. 
We also found that the suggested tetragonal deformation 
lifts the degeneracies. Thus we conclude that the cubic 
structure is unstable for negative J2. The lattice defor- 
mation may be spin-frustration-driven symmetry break- 
down like the spin Jahn- Teller (JT) effect 0. In MoaSbr, 
the spin frustration triggers the tetragonal deformation 
to decouple the 3D isotropic AF interaction, and stabi- 
lizes spatial long-range order of spin-singlet dimers, i.e., 
the valence bond crystal. As the reduced x from higher 
T Ts) and the considerably broadened specific heat 
(C) anomaly were reported 10], this state is realized in- 
stead of long-range AF order due to strong intra-dimer 
coupling already developed from sufficient high T. 

Taking account of the motional narrowing of nuclear 
dipolar fields observed by /iSR at ~ Ts, and structural 
characteristics of the low-T phase, we propose a sce- 
nario for the transition. First, note that the deformed 
states along [100], [010] and [001] directions are degen- 
erate. Here we refer to the corresponding Mo electronic 
states as ipx, ''Py and tpz- We expect that the high-T 
cubic state is the dynamic mixing of the three states, 
i.e., (V'k + ipy + ipz)/'^, just like the vibronic dynamic 
JT state. In other words, harmonic oscillations of the 
valence bonds keeping strong intra-dimer interaction are 
thermally induced. The narrowing of /xSR-A, which re- 
flects a relatively high-energy part of the dynamics, is 
explained naturally; at sufficiently low T, muons sense 
well defined nuclear dipolar fields from the deformed Mo 



sublattice, which are averaged by the dynamic atomic 
motion with increasing T. NQR-I/T2 above Ts, which 
nearly follows /iSR-A, is explained by the same mech- 
anism, while the damping below Ts is ascribed to the 
decoupling of the nuclear spin-spin interaction, i.e., the 
reduction of Sb like-nuclei due to the separation of the 
Sb sites in the distorted state. It is of interest that mag- 
netic entropy estimated from C exceeds i?ln2 ~ 0.69i? 
corresponding to the 5 = 1/2 manifold and is asymptotic 
to a much higher value ~ l.Oi? [10] close to i? In 3 ~ l.li? 
expected for the triply degenerate state. 

In conclusion, NQR and /iSR experiments revealed 
symmetry breakdown to a nonmagnetic state below Ts ~ 
50 K, with most probably tetragonal lattice distortion. 
Assuming negative Ji and J2 in the Mo sublattice, we 
pointed out possible frustration and proposed long-range 
order of spin-singlet dimers to the valence bond crystal, 
which transfers to the characteristic dynamic state with 
the cubic symmetry. Precise structural determination 
and experiments to detect phonon anomalies would be 
helpful to discuss the mechanism of the transition. On 
the other hand, NQR-l/Ti demonstrated the s-wave SC 
emerged from a metallic state where high-T Mo spins 
are quenched. The electron itinerancy and the present 
description focusing on the local nature should be recon- 
ciled in further studies, similarly to in other geometrically 
frustrated metals ^] . 
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